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Abstract 

We report all-electron and pseudopotential calculations of the ground-state energies of the neutral 
Ne atom and the Ne+ ion using the variational and diffusion quantum Monte Carlo (VMC and DMC) 
methods. We investigate different levels of Slater-Jastrow trial wave function: (i) using Hartrce-Fock 
orbitals, (ii) using orbitals optimized within a Monte Carlo procedure in the presence of a Jastrow factor, 
and (iii) including backflow correlations in the wave function. Small reductions in the total energy are 
obtained by optimizing the orbitals, while more significant reductions are obtained by incorporating 
backfiow correlations. We study the finite-time-step and fixed-node biases in the DMC energy, and 
show that there is a strong tendency for these errors to cancel when the first ionization potential (IP) 
is calculated. DMC gives highly accurate values for the IP of Ne at all the levels of trial wave function 
that we have considered. 
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1 Introduction 

Accurate approximations to the many-electron wave function are required as inputs for the variational and 
diffusion quantum Monte Carlo (VMC and DMC) electronic-structure methods The quality of these 
"trial" wave functions determines both the statistical efficiency of the methods and the final accuracy that 
can be obtained, and one of the main technical challenges in the field of quantum Monte Carlo (QMC) 
simulation is to develop more accurate trial wave functions. 

We report a detailed VMC and DMC study of the ground states of the Ne atom and the Nc+ ion, in 
which several different forms of trial wave function have been used, and both all-electron and pseudopotential 
calculations have been performed. The difference between the ground-state energies of the neutral atom and 
the positive ion gives the first ionization potential (IP), which is known accurately from experiments [2]. The 
IP's calculated using different wave functions arc compared with the essentially exact nonrclativistic IP [3] 
and the experimental value. We also demonstrate the degree to which the DMC fixed-node and time-step 
errors cancel when the IP is calculated. 

The basic form of wave function that we study consists of a product of Slater determinants for spin-up 
and spin-down electrons containing Hartree-Fock (HP) orbitals, multiplied by a positive Jastrow correlation 
factor. We consider the effects of optimizing the orbitals in the presence of the Jastrow factor [11 [5], and 
the effects of including backflow correlations in the wave function [SI [71 [H] . We have also carried out some 
tests using multidetcrminant wave functions [9l [10] , but we were not able to lower the DMC energy using 
this approach. 

The rest of this paper is organized as follows. In Sec. [H we briefly review the VMC and DMC methods, 
while in Sec. [H we describe the different levels of Slater- Jastrow wave function that we investigate. In Sec. [H 
we study time-step bias in DMC calculations of the total energy and IP of Ne. The total energies and IP's 
calculated using different trial wave functions are compared in Sec. [5l Finally, we draw our conclusions in 
Seel 

We use Hartree atomic units, ?i = |e| = rrie — 47reo = 1, throughout. All of our QMC calculations were 
performed using the CASINO code [lT| . 
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2 QMC methods 



The VMC energy is calculated as the expectation value of the Hamiltonian operator with respect to a 
trial wave function, the integrals being performed by a Monte Carlo method. In DMC, the imaginary- 
time Schrodinger equation is used to evolve a set of configurations towards the ground-state distribution. 
Fcrmionic symmetry is maintained by the fixed-node approximation |12| , in which the nodal surface of the 
wave function is constrained to equal that of a trial wave function. Our DMC algorithm is essentially that of 
Umrigar et al. jl3j . Various modifications to the DMC Green's function were proposed by Umrigar et al. to 
reduce the time-step errors in all-electron calculations, and we investigate the effects of these modifications 
in Sec. IH 

For most of our pseudopotential calculations we used a nonrelativistic HF Nc pseudopotential (effective 
core potential), although for our investigations of time-step errors we used a Dirac-Fock averaged- relativistic 
effective potential (AREP) [T31II1]. These pseudopotcntials arc finite at the origin, making them particularly 
suitable for use in QMC calculations. There is evidence [TB] that neglecting correlation effects entirely when 
constructing pseudopotcntials, as in HF and Dirac-Fock theory, provides better pseudopotcntials for use in 
correlated valence calculations than density-functional-theory-derived pseudopotcntials. Part of the motiva- 
tion for performing the calculations reported here is to test the Ne pseudopotcntials in correlated valence 
calculations. The contributions to the total energy from the nonlocal components of the pseudopotential 
within DMC were calculated using the locality approximation |17j . which leads to errors that arc second 
order in the error in the trial wave function |18j . 

3 Trial wave functions 

3.1 Basic Slater- Jastrow wave function with HF orbitals 

The Slater-Jastrow wave function may be written as 




(1) 
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where R = (ri, . . . ,rjv) denotes the spatial coordinates of all the electrons, exp[J(R)] is the Jastrow cor- 
relation factor, _dJ(R) and Dj'(R) arc Slater determinants of orbitals for spin-up and spin-down electrons, 
and the {fii} are expansion coefficients. The Jastrow factor, which describes the dynamic correlation of the 
electrons, is an explicit function of the interparticle distances. We use the form of Jastrow factor described in 
Ref. [in]. In our calculations the Slater determinants {Df} were formed from single-particle orbitals obtained 
from restricted open-shell HF calculations using numerical integration on radial grids. The all-electron HF 
calculations were performed using the code of Fischer [20j and the pseudopotential HF calculations were 
performed using our own code. The variable parameters in the trial wave functions were determined by 
minimizing the unrcweighted variance of the energy [H 1211 122j . 

We also performed tests using multideterminant trial wave functions. The determinants were obtained 
from multiconfiguration HF calculations performed by numerical integration on radial grids, the expansion 
coefficients [{fJ-i} in Eq. ([T])] being re-optimized in the presence of the Jastrow factor using unrcweighted 
variance minimization |22j. During the optimization process, we constrained the coefficients of the deter- 
minants within each configuration state function to have the same magnitude, so that the symmetry of the 
wave function was maintained. We were able to achieve significantly lower VMC energies than we could with 
the single-determinant Slater-Jastrow wave function, but the DMC energies were higher than the single- 
determinant ones. All numerical results reported in this article were obtained with a single term in the 
determinant expansion. 

3.2 Modification of the HF orbitals 

Altering parameters in the orbitals moves the nodal surface of the trial wave function and therefore changes 
the value of the fixed-node DMC energy. This is not the case for parameters in the Jastrow factor. Optimizing 
the orbitals in the presence of a Jastrow factor generally improves the nodal surface of the wave function 
and therefore lowers the DMC energy, which satisfies a variational principle with respect to errors in the 
nodal surface |23[ 124] . Unfortunately, optimizing parameters that affect the nodal surface by unrcweighted 
variance minimization is problematic. At each iteration, a minimization is performed using a fixed set of 
sampling points in the configuration space, and the nodal surface deforms continuously as the parameter 
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values are altered. In the course of this deformation, the nodal surface may pass through one of the fixed 
configurations, at which point the estimate of the unrcweighted variance of the local energy {'^^^H'ii, 
where H is the Hamiltonian operator) diverges. This makes it difficult to locate the true minimum of the 
unreweighted variance. It is possible to improve the situation by using the weight-limiting scheme described 
in Ref. |22| , which effectively removes configurations whose local energies are diverging from the set used to 
calculate the variance. Alternatively, the reweighted variance-minimization algorithm can be used, where 
each configuration is weighted by the ratio of the square of the current wave function to the probability 
density according to which the configurations were distributed, so that the calculated variance is an estimate 
of the actual variance of the energy. Reweighted variance minimization does not suffer from the problem of 
divergences in the variance but, nevertheless, the unreweighted algorithm gives significantly lower variational 
energies in practice, and we used this method for all the optimizations reported here. However, the difficulty 
in locating the global minimum of the variance with respect to parameters that affect the nodal surface 
makes it more difficult to optimize the orbitals for two different systems to the same level of accuracy. 
The atomic orbitals that we have used can be expressed in the form 



where we have used spherical polar coordinates, with r, 9, and (j) being the radial, polar, and azimuthal 
coordinates of the point r. The origin is chosen to lie at the nucleus. Y/m(0, </>) is the {I, m)th spherical har- 
monic, and Pnf {i") is the HF radial function. (Note that, in practice, we use appropriate linear combinations 
of spherical harmonics with a given I in order to construct real orbitals.) The modification to the HF radial 
function is written as 



where Ani, Bni and {c^nU ■ ■ ■ ,CNn{\ arc optimizablc parameters. Obviously we must have Ani > and 
Bni > 0, otherwise the orbital will be non-normalizable. The chosen form of the modification to the radial 
function has considerable variational freedom and decays exponentially at large distances. The HF orbitals 
minimize the energy in the absence of a Jastrow factor, although they arc no longer optimal when a Jastrow 
factor is included; however we expect Ap„;(r) to remain small in this case. After some experimentation we 
decided to use TV = 3 for the all-electron calculations and TV = 6 for the pseudopotential ones, as these 



^nim{v) = [p^f (r) + Ap„,(r)] r'Yi^{e,4>). 



(2) 




(3) 
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choices allowed us to achieve the lowest variational energies in our tests. 

Kato derived a condition on the spherical average of the wave function about a bare nucleus of atomic 
number Z at the origin |25| : 



where = |ri| and (^) is the spherical average of the many-body wave function about = 0. If the 
wave function is nonzero at the nucleus (as is the case for the ground-state wave functions studied in the 
present work), then it is possible to carry out accurate and stable QMC calculations if and only if the Kato 
cusp condition [Eq. is satisfied [T]. It is easy to show that a Slater- Jastrow wave function obeys the 
electron-nucleus Kato cusp condition if each orbital satisfies Eq. (|4]) individually and the Jastrow factor is 
cuspless at the nucleus. (It is preferable to impose the electron-nucleus cusp condition via the orbitals rather 
than the Jastrow factor [ISj-) The spherical averages of orbitals with I ^ are zero, and therefore these 
orbitals obey the Kato cusp condition automatically, but we must impose Eq. ([4]) on the radial parts of the 
s orbitals. This gives us the requirement that ci„o = —Zcono, where Z is the atomic number of the atom in 
an all-electron calculation and Z = for a calculation using a pseudopotcntial which is finite at the nucleus. 

3.3 Backflow correlations 

We have also investigated the effect of incorporating backflow correlations in the trial wave function. Classical 
backflow is related to the flow of a fluid around a large impurity, and its quantum analog was flrst discussed 
by Feynman [B] and Fcynman and Cohen [7] in the contexts of excitations in "^Ile and the effective mass 
of a ^Hc impurity in liquid "^Hc. Backflow correlations have previously been used in trial wave functions 
for fermion QMC simulations of two-dimensional [S] and three-dimensional |26] electron gases, and metallic 
hydrogen |27j . The original derivation of backflow correlations by Feynman was based on the idea of local 
current conservation, although they can also be derived from an imaginary-time evolution argument [8l I28j. 

The backflow correlations are described by replacing the electron coordinates {r;} in the Slater determi- 
nants of Eq. ([!]) by "quasiparticle" coordinates {x^}, defined to be 




(4) 




(5) 
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where F and G are functions that contain parameters to be determined by variance minimization, and r-y ~ 
\ri — rj|. In our work, F{rij,ri,rj) and G{rij ,ri,rj) consist of smoothly truncated polynomial expansions in 
ri, rj, and Vij, of the same form as the isotropic terms that occur in the Jastrow factor introduced in Ref. 

m- 

In our all-electron calculations the electron-nucleus cusp conditions are enforced via the orbitals in the 
Slater determinants, and in order to preserve the cusp conditions when backflow is included we ensure that F 
and G go smoothly to zero as an electron approaches the nucleus. We use pseudopotentials that are smooth 
at the nucleus and we therefore choose G to be smooth at the nucleus in our pseudopotential calculations. 
The electron-electron Kato cusp conditions are enforced by the Jastrow factor in our calculations, and the 
backflow functions are chosen to preserve these conditions. Our implementation of backflow correlations for 
systems containing many atoms will be described elsewhere |29) . 

4 Investigations of time-step errors 

We have studied the time-step bias in all-electron and Dirac-Fock AREP pseudopotential DMC calculations 
of the ground-state energies of the Ne atom and the Ne+ ion [Ml [15] . Figures [T] and [2] show the DMC 
energy obtained at different time steps, both with and without the modifications to the all-electron DMC 
Green's function proposed in Ref. [13]. The all-electron energies enter the small-time-step regime (where 
the bias in the DMC energy is nearly linear in the time step) at time steps less than about 0.005 a.u. The 
root-mean-square distance diffused at each time step of 0.005 a.u. is about 0.12 a.u., which is close to the 
Bohr radius of Ne {l/Z = 0.1 a.u.); this is the length scale relevant to the core (Is) electrons, suggesting 
that they are responsible for much of the time-step bias. The time-step bias remains small up to much larger 
time steps in the pseudopotential calculations, again implying that the Is electrons are responsible for most 
of the bias in the all-electron results. The contribution to the time-step bias from the core region is expected 
to be large, because the energy scale of the innermost electrons is greater than the energy scale of the outer 
electrons, although this is offset by the fact that the wave function is relatively accurate in the vicinity of 
the nucleus. 

It is clear from Figs. [T] and [D that the modifications to the all-electron DMC Green's function successfully 
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eliminate much of the time-step bias due to the innermost electrons at small time steps. On the other hand, 
the form of the time-step bias is almost the same for the Ne atom and the Ne+ ion, suggesting that the 
bias will cancel out when the IP is calculated, irrespective of whether the Green's function modifications are 
used. It can be seen in Fig. [3l which shows the IP as a function of time step, that this is indeed the case. 
Since the bias is largely caused by the innermost electrons, which are in almost identical environments in Ne 
and Ne+, this strong cancellation is to be expected. 

If a basic Slater- Jastrow wave function is used then the zero-time-step limit of the DMC energy is 
determined by the nodes of the HF wave function, independent of the Jastrow factor, whereas the DMC 
energy at finite time steps does depend on the Jastrow factor. However, provided that similar levels of 
Jastrow factor are used for systems to be compared (with the core electrons in the same configuration), it 
is reasonable to expect a large degree of cancellation of bias to occur when energy differences are taken at 
finite time steps. 

5 Energies of Ne and Ne+ 

In Table 1 we present values for the total nonrclativistic energies of Ne and Ne+ and the nonrelativistic 
IP of Ne, calculated using a number of different electronic-structure methods. For the all-electron atom 
and ion, we give results obtained using HF theory (AHF), density- functional theory using the local spin- 
density approximation (DFT-LSDA), Perdew-Wang (DFT-PW91), and Bccke-Lce- Yang-Parr (DFT-BLYP) 
exchange-correlation functionals, all-electron VMC and DMC using a basic Slater- Jastrow wave function 
with HF orbitals (AVMC and ADMC), all-electron VMC and DMC using a Slater- Jastrow wave function 
with optimized orbitals (AO VMC and AODMC), all-electron VMC and DMC with HF orbitals and backflow 
correlations (ABVMC and ABDMC), and all-electron VMC and DMC with optimized orbitals and backflow 
correlations (AOBVMC and AOBDMC). For the pseudopotential calculations, which used a nonrelativistic 
HF pseudopotential [Ml US], we give results obtained using HF theory (PHF), VMC and DMC using a basic 
Slater- Jastrow wave function with HF orbitals (PVMC and PDMC), VMC and DMC using a Slater- Jastrow 
wave function with optimized orbitals (POVMC and PODMC), VMC and DMC with HF orbitals and 
backflow correlations (PBVMC and PBDMC), and VMC and DMC with optimized orbitals and backflow 
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correlations (POBVMC and POBDMC). For each of the VMC calculations we give the total number of 
parameters in the trial wave function that were optimized using unreweighted variance minimization, along 
with the expectation value of the variance of the energy. Furthermore, we give the fraction of the correlation 
energy retrieved by the wave function. In the all-electron calculations, the correlation energy is defined to 
be the difference between the HF energy and the exact nonrelativistic energy. In the pseudo-Ne calculations, 
the correlation energy is defined to be the difference between the HF energy and the PBDMC energy. The 
DFT-LSDA IP is taken from Ref. 00], while the DFT-PW91 and DFT-BLYP resuhs are taken from Rcf. [31]. 
The "exact" nonrelativistic (NR) infinite- nuclear-mass energies are taken from Ref. [3] . 

All the DMC energies quoted in Table 1 have been extrapolated to zero time step. We used a range of 
small time steps and performed linear extrapolations of the energies to zero time step. 

Optimizing the orbitals reduces the VMC energy, but does not have a significant effect on the DMC 
energy, except in the case of the pseudopotential calculation with backfiow correlations, where optimizing 
the orbitals actually increases the DMC energy, presumably because the very large number of variable 
parameters adversely affects the optimization. The VMC variance is very slightly reduced by optimizing the 
orbitals in the all-electron calculations, but not in the pseudopotential calculations. (Note that, as discussed 
in Ref. |22| . the unreweighted variance-minimization algorithm does not generally minimize the expected 
variance of the energy.) The inclusion of backfiow correlations reduces the errors in the VMC energies by a 
factor of about 2, and reduces the variances of the VMC energies by a factor of between 2 and 3. For both 
pseudopotential and all-electron calculations, the percentage of the total correlation energy retrieved within 
VMC is slightly larger for Nc than for Ne"*", although the difference is smaller in DMC. 

The previous lowest VMC energy in the literature for the all-electron Ne atom is due to Huang et al. [52] , 
who used a single-determinant trial wave function with optimized orbitals and Jastrow factor, corresponding 
to our AOVMC level of calculation, and obtained an energy of —128.9029(3) a.u., which amounts to 91.11(8)% 
of the correlation energy. Our AOVMC energy of —128.90334(9) a.u. (91.23(2)% of the correlation energy) 
is very similar. We obtain a VMC energy of —128.9205(2) a.u. (95.62(5)% of the correlation energy) at the 
AOBVMC level, so that backfiow retrieves about 50.1(7)% of the remaining correlation energy. Note also 
that the VMC energies obtained with backfiow are very similar to the DMC energies without backfiow, for 
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both the aU-electron and pseudopotential calculations. This indicates that VMC with backflow may be a 
useful level of theory, because VMC calculations are significantly less costly than DMC ones, and VMC has 
advantages for calculating expectation values of quantities other than the energy. 

At the DMC level backflow retrieves about 38(2)% of the correlation energy missing at the AODMC level. 
Our AOBDMC energy of -128.9290(2) a.u. (97.80(5)% of the correlation energy) is significantly lower than 
the previous lowest DMC energy in the literature of —128.9243(8) a.u. (96.6(2)% of the correlation energy) 

m- 

The increase in the complexity of the wave function resulting from the inclusion of backflow correlations 
is apparent from the three- to five-fold increase in the number of parameters in the wave function. However, 
the use of backflow or orbital optimization does not result in any clear improvement to the VMC or DMC 
IP's. The all-electron DMC IP's obtained with the different forms of wave function are very close to one 
another. The inclusion of backflow in the pseudopotential calculations results in a small increase in the 
IP (away from the experimental result). This may reflect the fact that it is harder to optimize a backflow 
function for Ne^, where the symmetry between up and down spins is broken, or it may be an indication of 
the error inherent in the pseudopotential approximation. 

Irrespective of the form of trial wave function used, or whether all-electron or pseudopotential calculations 
are performed, or whether the VMC or DMC method is used, all of the QMC calculations give excellent 
values for the IP. In the worst cases the errors in the VMC nonrelativistic IP's are 0.43(4)% (all-electron) 
and 0.60(3)% (pseudopotential) and the errors in the DMC nonrelativistic IP's are 0.28(4)% (all-electron) 
and 0.35(4)% (pseudopotential). The considerably larger errors in the HF IP (8%) and DFT IP's (1.7% 
for the best case of the BLYP functional) demonstrate the importance of describing electron correlation 
accurately when calculating the IP. 

The correction to the IP of Ne from the finite mass of the nucleus is of order 10^^ a.u., which is negligible 
on the scale of interest. The total relativistic correction to the IP of Ne has been estimated to be —0.00196 a.u. 
[3], which is significant on the scale of interest: see Fig. [31 Adding this correction to our best DMC data, 
obtained at the AOBDMC level, gives an IP of 0.7945(2) a.u., which is only 0.0020(2) a.u. (0.054(5) eV) 
larger than the experimental value of 0.792481 a.u. (21.5645 eV) [2]. Using the same relativistic correction for 
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the pseudopotential POBDMC data gives an IP which is 0.0026(3) a.u. (0.071(8) eV) larger than experiment. 
The error in the IP due to the use of the pseudopotential is therefore about 0.0006(4) a.u. (0.016(11) eV), 
which is less than 0.1% of the IP. 

6 Conclusions 

We have performed all-electron and pseudopotential VMC and DMC calculations of the total energies of 
Ne and Ne+ using basic Slater- Jastrow wave functions with HF orbitals, Slater- Jastrow wave functions with 
optimized orbitals, and wave functions including backflow correlations. We have found that optimizing the 
orbitals makes a small improvement to the VMC energy, but has very little effect on the DMC energy. The HF 
orbitals give nearly optimal single-determinant nodal surfaces for the ground states of Ne and Ne+, although 
this cannot be expected to hold for other systems. On the other hand, including backflow correlations lowers 
both the VMC and DMC energies significantly, giving us the lowest VMC and DMC energies for all-electron 
Ne in the literature to date. However, the improvements in the total energies of Ne and Ne"*" hardly affect 
the calculated IP, because they almost exactly cancel when the energy difference is evaluated. Overall, the 
calculated IP's are slightly too large, because the wave functions for Ne+ are not quite as accurate as those 
for Ne. 

Including backflow correlations for Ne and Ne+ reduced the error in the VMC total energy by a factor of 
about 2 and reduced the variance of the VMC energy by a factor of between 2 and 3. Reducing the variance 
improves the statistical efficiency of VMC calculations, while reducing the error in the VMC total energy 
improves the statistical efficiency of DMC calculations [33l [34] . The incorporation of backflow is expected to 
improve the QMC estimates of all expectation values, not just the energy. Including backflow correlations 
is costly, however, because every element in the Slater determinant has to be recomputed each time an 
electron is moved, whereas only a single column of the Slater determinant has to be updated after each 
move when the basic Slater-Jastrow wave function is used. In the case of Ne, including backflow correlations 
increased the computational cost per move in VMC and DMC by a factor of between 4 and 7. However, 
the reduction in the variance meant that the number of moves required to obtain a fixed error bar in the 
energy was smaller; hence, overall, including backflow correlations increased the time taken to perform the 
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calculations by a factor of between 2 and 3. Backflow functions introduce additional variable parameters into 
the trial wave function, making the optimization procedure more difficult and costly. Nodal surfaces obtained 
from HF or DFT orbitals do not suffer from statistical errors. By contrast, the parameters in the backflow 
functions are optimized by a Monte Carlo method, and therefore the nodal surface is subject to stochastic 
noise. Furthermore, although Jastrow factors are easy to optimize using unreweighted variance minimization 
with a fixed distribution of configurations, parameters that affect the nodal surface are relatively difficult to 
optimize, because the local energy diverges at the nodes. 

We have investigated the time-step bias in the DMC total energies of Ne and Ne+ and the IP of Ne. The 
time-step errors in the pseudopotcntial calculations were relatively small, indicating that the bias in the all- 
electron total energy was mainly due to the Is electrons. In IP calculations using a basic Slater- Jastrow wave 
function with HF orbitals the time-step bias largely canceled out when the difference of energies was taken, 
irrespective of whether the modifications to the all-electron DMC Green's function proposed in Ref. [13] were 
used. The time-step errors in the IP were considerably smaller in the pseudopotential calculations than the 
all-electron ones. 

Most importantly, our results demonstrate the superb accuracy of the DMC method for this problem. Af- 
ter including a correction for relativistic effects, our all-electron IP of Ne is about 0.0020 a.u. (0.054 eV) larger 
than the experimental value [2] of 0.792481 a.u. (21.5645 eV). Our pseudopotential IP is about 0.0026 a.u. 
(0.071 eV) larger than the experimental value, demonstrating the accuracy of the pseudopotentials we use 

[Hill]- 
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Table 





No. 


par am. 


Energy 


(a.u.) 


%age 


corr. en. 


Variance (a.u.) 




Meth. 


Ne 


Ne+ 


Nc 


Ne+ 


Ne 


Ne+ 


Nc 


Ne+ 


IP (a.u.) 


Exact NR 




- 


-128.9376 


-128.1431 


100.0 


100.0 


- 


- 


0.7945 


DFT-LSDA 


- 


- 


- 


- 


- 


- 


- 


- 


0.812 


DFT-PW91 


- 


- 


- 


- 


- 


- 


- 


- 


0.812 


DFT-BLYP 


- 


- 


- 


- 


- 


- 


- 


- 


0.808 


AHF 


- 


- 


-128.54710 


-127.81781 


0.0 


0.0 


- 


- 


0.72928 


AVMC 


75 


79 


-128.8983(2) 


-128.1016(1) 


89.94(5) 


87.24(3) 


1.12(2) 


1.047(1) 


0.7967(2) 


AOVMC 


88 


92 


-128.90334(9) 


-128.10760(6) 


91.23(2) 


89.09(2) 


0.9740(9) 


0.938(3) 


0.79574(11) 


ABVMC 


205 


209 


-128.9190(2) 


-128.1213(2) 


95.24(5) 


93.30(6) 


0.421(1) 


0.406(2) 


0.7979(3) 


AOBVMC 


218 


222 


-128.9205(2) 


-128.1240(2) 


95.62(5) 


94.13(6) 


0.391(1) 


0.379(3) 


0.7965(3) 


ADMC 


- 


- 


-128.9233(2) 


-128.1266(2) 


96.34(5) 


94.93(6) 


- 


- 


0.7967(3) 


AODMC 


- 


- 


-128.9237(2) 


-128.1273(2) 


96.44(5) 


95.14(6) 


- 


- 


0.7964(3) 


ABDMC 


- 


- 


-128.9287(2) 


-128.1321(1) 


97.72(5) 


96.62(3) 


- 


- 


0.7966(2) 


AOBDMC 


- 


- 


-128.9290(2) 


-128.1325(1) 


97.80(5) 


96.74(3) 


- 


- 


0.7965(2) 


PHF 


- 


- 


-34.59111 


-33.86129 


0.0 


0.0 


- 


- 


0.72982 


PVMC 


53 


61 


-34.88930(4) 


-34.09453(4) 


93.10(6) 


92.26(7) 


0.3955(6) 


0.3994(8) 


0.79477(6) 


POVMC 


70 


76 


-34.88973(4) 


-34.09472(4) 


93.23(6) 


92.33(7) 


0.4024(7) 


0.391(1) 


0.79501(6) 


PBVMC 


275 


291 


-34.90509(6) 


-34.1059(2) 


98.03(6) 


96.8(1) 


0.1748(7) 


0.226(2) 


0.7992(2) 


POBVMC 


292 


300 


-34.90475(6) 


-34.1055(2) 


97.92(6) 


96.6(1) 


0.1483(7) 


0.222(2) 


0.7993(2) 


PDMC 






-34.9026(2) 


-34.1072(2) 


97.25(9) 


97.3(1) 






0.7954(3) 


PODMC 






-34.9028(2) 


-34.1074(2) 


97.31(9) 


97.3(1) 






0.7954(3) 


PBDMC 






-34.9114(2) 


-34.1141(2) 


100.00(9) 


100.0(1) 






0.7973(3) 


POBDMC 






-34.9099(2) 


-34.1128(2) 


99.53(9) 


99.5(1) 






0.7971(3) 



Tabic 1: Numbers of parameters in trial wave functions, ground-state energies, percentages of correlation 
energies retrieved, energy variances, and IP's of all-electron Nc (and Nc+) and pseudo-Ne (and pseudo-Ne+), 
obtained using various methods. 
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Figure captions 

1. DMC energy of the all-electron Ne atom and the pseudo-Ne atom as a function of the time step. The 
pseudo-Ne results are offset so that the zero-time-step energy matches the all-electron energy. The 
statistical error bars are smaller than the symbols. "GF mods" refers to the modifications to the DMC 
Green's function proposed in Ref. [13j . The "exact" nonrelativistic infinite nuclear mass result is taken 
from Ref. [3]. 

2. The same as Fig. [1] but for a Ne"*" ion. 

3. IP of all-electron and pseudo Ne as a function of the time step. "GF mods" refers to the modifications 
to the DMC Green's function proposed in Ref. [13]. The experimental result is taken from Ref. [2], and 
the exact nonrelativistic IP is from Ref. [3] . The statistical error bars are smaller than the symbols. 
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Figure [T] 
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Figure [2] 
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Figure [3] 
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